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ABSTRACT 


This  paper  describes  structural  reliability  developed  from  an  assumed  flaw 
simulation  scheme  combined  with  an  idealized  linear  elastic  fracture  mechanics 
model.  In  the  reliability  computation,  it  has  been  demonstrated  that  application 
of  this  scheme  can  provide  a  sensitivity  analysis  in  relation  to  flaw  detection 
capabilities.  To  demonstrate  the  versatility  of  the  method,  reliability  numbers 
were  obtained  for  both  a  simply  loaded  fragmentation  shell  and  an  antitank  projec¬ 
tile  subjected  to  a  complex  stress  state.  The  flaw  simulation  scheme  results  are 
compared  to  corresponding  reliability  determinations  from  conventional  Weibull 
and  Warner  stress-strength  diagram  methods. 

The  reliability  definitions  associated  with  the  scheme  require  a  probability 
density  function  representation  of  the  material  strength  and  stress  distributions 
for  each  structural  element.  A  density  function  representation  of  the  allowable 
stress  (strength)  is  obtained  from  fracture  mechanics  Kjq  relationships  in  con¬ 
junction  with  the  Monte  Carlo  method  where  a  specific  random  form  of  the  parame¬ 
ter  is  assigned.  Crack  orientations  are  assumed  to  vary  in  a  uniform  random 
manner  with  respect  to  the  principal  axis  and  its  normal.  The  crack  size  varia¬ 
tion  is  defined  in  exponential  functional  form  where  the  sizes  vary  from  a  large 
number  of  very  small  cracks  to  a  relatively  small  number  of  larger  cracks.  The 
structural  configuration  determines  the  type  of  crack  and  its  location.  It  is 
assumed  that  flaws  exist  in  the  structural  elements  and  are  remote  from  any  other 
flaws.  A  normal  density  function  represented  element  stress  variability  where 
the  mean  stress  was  obtained  from  an  axisymmetric  analysis  using  a  unique  finite 
element  code.  Perturbations  of  the  coefficient  of  variation  for  the  assumed  nor¬ 
mal  distribution  monitored  the  effects  of  errors  in  the  finite  element  solution. 

Both  the  "weakest  link"  and  series-parallel  system  are  evaluated  for  desir¬ 
ability  in  estimating  structural  reliability.  The  weakest  link  approach,  which 
introduces  reliability  independence  between  elements,  will  thereby  describe  a 
much  more  conservative  reliability  estimate  than  the  series-parallel  system  which 
requires  at  least  two  adjacent  elements  to  fail  in  order  to  have  structural 
fai  lure . 

Examining  the  reliability  computation  for  both  structures  as  a  function  of 
assumed  minimum  detectable  crack  sizes  indicated  the  importance  of  determining 
relatively  small  cracks  to  obtain  acceptable  reliability. 
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INTRODUCTION 


Present  linear  elastic  fracture  mechanics  (LEFM)  procedures  as  applied  to  structural 
reliability  depend  on  determination  of  critical  crack  sizes  by  using  proper  Kjq  (fracture 
toughness)  relationships  in  conjunction  with  known  stress  distributions  in  the  structure. I 
In  obtaining  critical  size,  it  is  assumed  the  cracks  are  orientated  normal  to  the  maximum 
principal  stress.  Once  the  critical  crack  sizes  and  locations  are  established,  a  non¬ 
destructive  test  (NDT)  procedure  is  applied  to  the  structure  to  determine  if  this  crack 
exists.  If  a  critical  crack  is  located,  the  structure  is  rejected.  It  should  be  noted 
that  an  excessively  large  rejection  rate  can  occur  by  applying  this  method  since  the 
probability  of  the  assumed  crack  orientation  is  very  small. 2  Another  objection  to  the 
LEFM  procedure  is  the  inability  of  NDT  methods  to  detect  cracks  less  than  a  certain 
size  (e.g.,  0.10  inch)."  If  the  critical  size  is  computed  to  be  less  than  this  size 
then  the  LEFM  procedure  may  not  be  able  to  determine  potential  failure  in  the  structure. 

The  uncertainties  existing  in  the  use  of  the  analytic  tools,  such  as  finite  element 
(FE)  analysis,  in  obtaining  the  structural  stress  distribution  can  introduce  consider¬ 
able  error  in  obtaining  the  critical  flaw  size.  For  example,  the  accuracy  of  the  FE 
solution  depends  on  the  severity  of  stress  gradients,  mesh  size,  types  of  elements 
used,  and  the  effects  of  averaging  stresses  within  the  element.  Therefore,  it  is  impor¬ 
tant  to  recognize  element  stress  variability  in  establishing  structural  reliability. 

The  flaw  simulation  scheme  (FSS)  introduced  in  this  paper  attempts  to  provide 
alternatives  to  the  conventional  methods  described  above.  The  FSS  procedure  does  not 
describe  exact  reliability  numbers  but  rather  determines  the  effects  of  flaw  detection 
capabilities  and  stress  computation  errors.  Applying  this  sensitivity  analysis  also 
determines  which  possible  parametric  errors  most  affect  the  reliability  computations. 

For  example,  if  Kjq  numbers  obtained  from  laboratory  test  results  indicated  a  15  per¬ 
cent  coefficient  of  variation  (CV)  and  the  resultant  reliability  computation  was  0.850 
versus  0.999  for  an  assumed  CV  of  5  percent,  then  one  should  attempt  to  reduce  the 
variability  in  this  parameter. 

In  using  the  Monte  Carlo  method, ^  the  flaw  simulation  scheme  provides  for  variations 
in  crack  orientation  and  size  in  addition  to  computed  stress  values  in  the  structure. 
Crack  orientations  are  assumed  to  vary  from  0°  to  90°  in  a  uniform  random  manner.  The 
crack  size  variation  is  defined  in  exponential  functional  form  where  size  varies  from  a 
large  number  of  very  small  cracks  to  a  relatively  small  number  of  larger  detectable 
cracks.^  The  types  of  cracks  and  their  locations  are  the  through-center  crack,  near 
cut-out  edge,  corner  crack,  and  surface  crack  (center).  The  structural  configuration 
determines  the  types  of  cracks  and  their  locations. 

The  stress  values  obtained  for  a  cracked  structural  element  is  assumed  to  be  a 
normal  distribution  where  the  CV  is  varied  in  order  to  determine  the  effect  of  errors 
in  the  FE  analysis. 
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The  FE  code  used  in  obtaining  stress  values  utilizes  an  8-noded  isoparametric 
element  with  substructuring  capabilities.  Each  element  contains  16  degrees  of  freedom 
with  the  complete  structure  represented  by  a  minimum  of  3000  degrees  of  freedom.  In 
the  substructuring  procedure,  the  stiffness  matrix  solutions  are  modified  in  the  reduc¬ 
tion  process  to  provide  selected  independent  determination  of  the  displacement  coeffi¬ 
cients  . 

The  loading  conditions  applied  to  the  structures  used  in  demonstrating  the  applica 
bility  of  the  FSS  are  only  approximations.  Therefore,  one  should  not  assume  any  level 
of  quality  assurance  of  either  the  fragmentation  shell  or  the  antitank  projectile  (ATP) 

TYPES  OF  STRUCTURES  CONSIDERED 

The  model  problems  chosen  to  illustrate  the  technique  proposed  in  this  paper  are 
the  determination  of  structure  reliability  for  the  ATP  control  section  (Figure  1)  and  a 
fragmentation  shell  (Figure  2).  Both  structures  are  made  from  relatively  high  strength 
brittle  steels.  The  projectile  control  section  is  subjected  to  relatively  large  com¬ 
pressive  loads  at  the  aft  section  and  tensile  stresses  in  the  vicinity  of  the  cut-out 
region  where  the  fins  are  attached.  The  fragmentation  shell  is  internally  loaded  with 
a  uniform  pressure  of  14  ksi.  This  load  represents  a  proof  test  used  in  evaluating 
shell  quality.  Reliability  determination  of  the  control  section  provided  the  opportu¬ 
nity  to  evaluate  the  series-parallel  system  approach  because  of  the  complex  stress 
state  introduced  by  the  relatively  large  tension  and  compressive  stresses  in  the 
structure.  The  weakest  link  approach  is  more  readily  adapted  to  the  fragmentation 
shell  loading  state. 


Figure  1.  Antitank  projectile  configuration. 


Figure  2.  Cross-sectional  view  of 
fragmentation  shell.  (1)  Ogive, 
(2)  Bourrelet,  (3)  Inside  Radius, 
(4)  Thread  Relief,  and  (5)  Rear 
of  Base. 
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STATISTIC  EVALUATION  OF  VARIABLE  STRENGTH 


The  lack  of  ductility  characteristic  of  brittle  materials  has  two  undesirable 
consequences  for  the  engineer.  First,  any  misfit  or  misalignment  produces  locally  high 
stresses  which  cannot  be  relieved  by  plastic  flow,  unlike  ductile  materials.  Brittle 
component  designs  differ  from  those  for  similar  ductile  components  in  that  extra  atten¬ 
tion  must  be  paid  to  detail,  especially  in  highly  stressed  areas.  The  second  conse¬ 
quence  is  more  fundamental;  all  materials  contain  flaws  such  as  microscopic  cavities 
and  dislocations.  In  loaded  brittle  materials,  these  flaws  result  in  local  stress  con¬ 
centration  within  the  material.  The  strength  of  a  component  is  governed  by  the  chance 
that  a  severe  stress  concentration  (c)  will  be  subjected  to  a  stress  (c)  such  that  the 
local  stress  co  exceeds  the  material  strength.  The  occurrence  of  this  is  a  matter  of 
chance  and  explains  the  marked  variability  generally  observed  in  brittle  material 
strengths.  It  also  explains  why  brittle  material  failures  may  start  away  from  the 
maximum  stress,  and  failure  may  occur  at  a  severe  flaw  subject  to  a  lower  stress  at  a 
position  where  co  is  a  maximum. 

To  overcome  the  strength  variability  by  drastically  reducing  the  applied  loading 
is  not  an  attractive  engineering  proposition  What  is  needed  is  an  estimate  of  the 
likelihood  of  failure  of  the  component  under  a  specified  load.  This  requires  a  de¬ 
tailed  knowledge  of  the  stresses  in  the  structure  and  the  flaw  distribution  in  the 
material.  Well-established  techniques  are  available  for  the  stress  analysis,  some  of 
which  are  mentioned  later.  The  variation  in  material  strength  due  to  flaws  in  a  par¬ 
ticular  material  can  be  illustrated  by  fracture  tests  on  a  sample  of  specimens.  A 
histogram  of  the  fracture  stresses  of  both  brittle  and  ductile  material  subjected  to 
uniform  tension  is  shown  in  Figure  3a;  the  frequency  of  failure  (Ff)  is  the  fraction  of 
the  sample  failing  within  the  stress  range  a  to  a+6cs.  In  the  limit,  as  the  number  of 
specimens  (N)  becomes  large,  the  stress  interval  (<$o)  in  Figure  3b  can  be  reduced  to 
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give  a  continuous  distribution  curve.  Note  the  relatively  barge  variation  in  strength 
of  the  brittle  material  as  compared  to  the  corresponding  ductile  material  of  similar 
test  specimens.  Structures  with  large  variations  in  material  strength,  as  shown  for 
the  brittle  material  strength,  require  a  probabilistic  approach  in  the  design  pro¬ 
cedures. 

A  complementary  form  of  Figure  3a  is  obtained  if  the  data  is  presented  in  terms  of 
the  cumulative  failure  probability  (Pf).  This  quantity  is  the  fraction  of  the  sample 
failing  at  or  below  the  stress  o;  in  the  limit  it  is  the  integral  of  the  frequency  dis¬ 
tribution  with  respect  to  stress,  i.e., 

PfOO  -/%  do.  (1) 

o 

In  practice,  the  estimated  cumulative  failure  probability  is  usually  found  from  the 
data  using  the  "mean  ranking"  approach.  The  N  failure  stresses  of  the  sample  are 
arranged  in  ascending  order:  the  cumulative  failure  probability  associated  with  the 
icfl  failure  stress  in  the  list  is 

Pf(a.)  =  i/N+1 .  (2) 

The  probability  distribution  of  the  data  can  be  plotted  from  this  (see  Figure  3b). 

CONVENTIONAL  FRACTURE  MECHANICS  APPROACH  TO  STRUCTURAL  RELIABILITY 

Fracture  mechanics  in  the  design  process  requires  the  consideration  of  three  factor: 
a  stress  analysis,  a  measure  of  Kjq,  and  the  capability  of  inspecting  for  cracks.  The 
stress  analysis  can  require  elaborate  analysis  using  advanced  FE  or  simpler  closed  form 
solutions  depending  on  loading  conditions  and  structural  geometry.  Determinations  of 
the  plane  strain  fracture  toughness  of  a  material  are  necessary  under  a  sufficiently 
wide  variety  of  conditions  to  allow  realistic  assessment  of  the  minimum  range  of  values 
likely  to  be  encountered  in  design  conditions.  Some  of  the  crack  detection  techniques 
include  ultrasonics,  dye  penetrants,  magnetic  particles,  and  visual  inspection. 

Analytically,  the  critical  flaw  size  is  defined  by  an  expression  of  the  type  shown 
below: 

ac  =  f  (Q’  KIC>  a)  (3) 

Where:  ac  =  critical  flaw  depth, 

Q  =  a  parameter  which  takes  into  account  the  shape  of  the  flaw, 

Kjc  =  plane  strain  fracture  toughness  of  the  material,  and 

o  =  the  tensile  component  of  stress  acting  normal  to  the  plane 
of  the  flaw. 

The  critical  flaw  depth,  on  the  basis  of  the  fracture  toughness  and  stress  factors, 
can  result  in  catastrophic  crack  propagation.  It  should  be  noted  that  o,  as  defined 
above,  assumes  the  flaw  is  normal  to  the  acting  stress  (see  Figure  4).  This  assumption 
rejects  the  possibility  that  the  flaw  could  be  oriented  in  other  directions,  thereby 
neglecting  obvious  possibilities  in  favor  of  an  unlikely  one.  This  could  result  in  in¬ 
correctly  determining  critical  flaw  size  ac .  Present  flaw  detection  methods  in  many 
instances  are  not  capable  of  detecting  critical  flaws  of  a  relatively  small  size.  In 
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Figure  4.  Conventional  single-valued 
crack  orientation. 


fact,  comments  by  Hastings  indicate  flaws  less  than  0.1  inch  may  not  be  found  in  a 
structure.  Of  course,  the  other  difficulties  involve  not  being  able  to  find  a  detect¬ 
able  flaw  in  a  structure  although  it  does  exist. 


FLAW  SIMULATION  METHOD 

An  alternative  to  the  previously  described  conventional  fracture  mechanics  approach 
is  made  by  introducing  variation  in  crack  orientation  and  length  by  means  of  the  Monte 
Carlo  method.  Initially,  four  types  of  cracks  are  to  be  considered  in  an  element  (see 
Figure  5).  The  type  and  location  of  cracks  depends  on  the  structural  configuration. 

In  the  simulation  scheme  the  allowable  uniform  stress  a  from  Kjq  relationships  is 
written  as: 

°c  =  f(KIC>  A»  6)  (4) 

where:  2  =  crack  length  and 

9  =  angle  of  inclination  of  crack  (see  Figure  6d). 

It  is  assumed  that  crc  represents  the  material  strength  and  depends  on  the  parameters 
KlC>  2)9-  The  variations  in  Kjq  are  represented  by  a  normal  probability  density  function 
(PDF)  (Figure  6a).  The  angle  9  is  represented  by  uniform  random  numbers  in  the  range 
of  0°  to  90°  (Figure  6b).  The  distribution  of  sizes  9  is  of  an  exponential  PDF  form 
shown  in  Figure  6c.  The  oc  distribution  is  obtained  from  generating  a  set  of  uniform 
numbers  and  solving  for  x  in  the  relation, 

/Xf  i  =  R’  (5) 

—  OO 

where  R  =  uniform  random  number  and  fj^  corresponds  to  the  desired  type  of  frequency 
distribution.  The  Kjq  distribution  requires  test  results  for  the  material  used  in  the 
structure  in  order  to  obtain  the  necessary  mean  and  standard  deviation  values  In 
Figure  6c,  the  maximum  crack  length  £2  is  represented  by  the  smallest  detectable  crack 
consistent  with  the  capability  of  presently  available  NDT  methods.  The  assumed  distri¬ 
bution  of  crack  sizes  has  been  substantiated  by  McClintock4  and  Smith." 

*SM!TH,  J.  M.  Army  Materials  and  Mechanics  Research  Center,  Watertown,  Massachusetts,  NDT  Industrial  Applications  Branch, 

Personal  Communication. 

4.  McCLINTOCK,  F.  A.  Statistics  of  Brittle  Fracture.  Fracture  Mechanics  of  Ceramics,  v.  1,  Plenum  Press,  New  York,  1973,  p.  93-114. 
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Figure  5.  Crack  types  assumed  in  simulation  scheme. 
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Figure  6.  Simulation  method. 


A  PDF  can  be  obtained  for  the  allowable  stress  ctc  by  randomly  selecting  from  Kjq, 
9  and  £  distributions  discrete  sets  of  numbers  and  substituting  them  into  Equation  4. 
Note,  there  should  be  an  equal  amount  of  N  random  numbers  for  each  parameter  to  have  N 
numbers  representing  the  oc  distribution. 
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The  Kjc  relationships  for  ac  are  written  as: 


(a)  Through  crack 

a  =  K  (2/tOi/20  sin^)-1/2 

C  J.  C 

where:  KjC  =  (K2  +  Kjj)1^2,  (6) 

Kj  =  a  (ilTr/2)  1/2sin20 ,  and 
Kjj  =  a(2tr/2)  i/2sin0cos0. 

(b)  Corner  crack 

oc  =  KIC/ (2(1.28) (Z/n)l/2).  (7) 

(c)  Surface  crack 

ac  =  KIC/l.l(7ra/Q)l/2 

where:  Q  =  E(K)2  -  0.212 (a/a  )2, 

ys 

E(K)  =  complete  elliptical  integral  of  the  second  kind, 

K2  =  1  -  (a/2.)2,  (g) 

0  <  a  <  0.10  inch,  and 
0  <  2  <  0.05  inch. 


The  oc  relationship  for  the  inclined  edge  crack  was  obtained  from  application  of  a 
modified  mapping  collection  scheme. ^  The  results  are  presented  in  Table  1.  The  appro¬ 
priate  interpolation  procedure  was  applied  to  use  the  Monte  Carlo  method  as  outlined  in 
Equation  4. 


Table  1.  INCLINED  EDGE  CRACK  SOLUTION 


e 

Hr  =  Kj/AT 

Hq  =  Kjj//jrT 

0 

0.000 

0.000 

10 

0. 160 

0.170 

20 

0.296 

0.280 

30 

0.461 

0.335 

40 

0.680 

0.372 

45 

0.705 

0.365 

50 

0.781 

0.354 

60 

0.920 

0.305 

70 

1.028 

0.224 

80 

1.098 

0.118 

90 

1.124 

0.000 

°c  = 

KIC/Hn  ^ 

*FREESE,  C.  E.  Army  Materials  and  Mechanics  Research  Center,  Watertown,  Massachusetts,  unpublished  results. 
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STRUCTURAL  STRESS  ANALYSIS 


A  finite  element  analysis  was  applied  to  obtain  the  stress  distributions  in  the 
ATP  and  fragmentation  shell.  The  loading  consists  of  a  set-back  compressive  load  acting 
at  the  base  of  the  ATP  and  an  internal  pressure  (proof-test  load)  applied  to  the  shell. 
Rectangular  elements  were  used  in  the  analysis  for  both  structures  where  the  shell 
contains  693  elements  and  the  ATP  contains  601  elements.  The  FE  solution  determines 
the  average  maximum  and  minimum  principal  stresses  in  each  element.  The  maximum  stress 
is  used  in  the  reliability  determinations.  These  stresses  should  not  be  confused  with 
critical  stresses  obtained  from  the  Kjq  relationships  previously  described.  Having 
obtained  the  element  stresses  and  corresponding  crc  or  strength  values,  reliability  of 
the  elements  can  be  determined. 


ELEMENT  RELIABILITY  CALCULATION  METHODS 

Element  reliability  as  related  to  the  stress-strength  (Warner)  diagram-*  method 
(see  Appendix  A)  assumes  that  the  probability  of  survival  (reliability)  is  the  probabil¬ 
ity  that  material  strength  will  be  greater  than  the  stress  in  a  given  structural  element 
over  a  range  of  stress  values.  The  uncertainties  in  the  FE  solution  are  represented  by 
a  normal  distribution  f2,  where  the  calculated  mean  stress  and  the  assumed  CV  are  the 
functional  parameters.  The  distribution  f]  is  obtained  from  known  strength  data  (e.g., 
laboratory  tests).  This  distribution  does  not  necessarily  have  to  be  a  normal  function. 
It  can  be  any  PDF  that  accurately  represents  the  empirical  ranked  data.  In  Appendix  A, 
the  probability  of  S|  occurring  is  f 2 ( S  | )  ,  while  the  probability  of  strength  greater 
than  S ^  is  represented  by  the  integral,  with  fj  the  integrand  and  limits  of  Sj  and 
Multiplication  of  these  elements  provides  the  necessary  independence  between  the  two 
conditions.  Finally,  integration  over  the  entire  range  of  stress  values  defines  proba¬ 
bility  of  survival  Ps  of  each  element. 

Element  reliability  numbers  were  also  obtained  from  an  approach  similar  to  the  one 
previously  described  except  that  discrete  values  obtained  from  FSS  were  used  to  repre¬ 
sent  oc .  Element  design  stress  is  normally  represented  in  the  same  manner  except  dis¬ 
tributed  numbers  are  obtained  from  the  new  stress  and  assumed  CV.  This  is  a  reasonable 
approach  since  distribution  of  strength  values  do  not  necessarily  conform  to  any  known 
density  function.  The  tails  of  the  two  density  functions  are  more  accurately  repre¬ 
sented  than  by  some  crude  approximating  function.  This  method  is  outlined  in  Appendix 
B  where  the  probability  of  element  survival  is  defined  as  follows:  =  1 ,  when 

strength  is  greater  than  stress  values,  otherwise  it  is  zero.  This  process  is  completed 
when  all  combinations  are  considered.  The  relationship  Psj^  therefore  defines  the  Kth 
element  reliability  number. 

The  Weibull  statistic  was  used  for  comparison  purposes  in  obtaining  element  reli¬ 
ability  of  the  fragmentation  shell.  This  procedure  is  recommended  for  brittle  materials 
subjected  to  the  uniform  tension  state,  which  is  consistent  with  the  structural  loads 
and  material  described  for  the  shell.  The  Weibull0  PDF  is  commonly  used  in  ceramics 
and  other  brittle  materials  evaluation.  It  employs  the  weakest  link  concept  which  is 
consistent  with  failure  phenomenon  of  brittle  materials  subjected  to  tensile  stresses 
primarily.  A  plot  of  strength  versus  cumulative  density  function  (CDF)  for  HF-I  steel 
used  in  the  fragmentation  shell  construction  is  shown  in  Figure  7.  Note  the  excellent 
correlation  between  empirical  data  and  the  Weibull  function. 

5.  HAUGEN,  E.  B.  Probabilistic  Approaches  to  Design.  John  Wiley  and  Sons,  Inc.,  New  York,  1968. 

6.  WEIBULL,  W.  A  Static  Theory  of  the  Strength  of  Materials.  Proceedings  of  the  Royal  Swedish  Academy  of  Engineering  Sciences, 
v.  151,  January  1939. 
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Stress  (ksi) 


Figure  7.  Cumulative  density  function 
versus  stress. 


The  Weibull  probability  of  survival  Pg .  for  individual  stressed  components  is 
written  as 


P  =  exp 

i 

where:  K  =  1  for  simple  tensile  stress, 

V£  =  volume  of  element, 

V*  =  volume  of  test  specimen, 
amaxi  =  maximum  principal  stress  in  the  element,  and 

ct0  and  m  are  functional  parameters  obtained  from  test  data  using  the  maximum  likelihood 
method.  It  is  obvious  that  Pg .  is  volume  dependent,  that  is,  larger  volume  provides 
smaller  P„.  numbers. 


STRUCTURAL  RELIABILITY 

In  order  to  obtain  Ps  of  the  entire  structure,  the  weakest  link  concept  is  applied, 
that  is,  it  is  assumed  that  each  event  or  element  probability  of  survival  is  independent 
of  any  other  one  in  the  structure.  Therefore,  the  total  Ps  is  written  as 

N 

p  =  n  Psv,  do) 

S  K=1  K 

where:  N  =  number  of  elements  in  the  structure  and 

PsK  =  the  probability  of  survival  of  the  individual  element. 

The  corresponding  probability  of  failure  is  defined  as  Pf  =  1-PS. 
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A  series-parallel  approach  is  introduced  to  examine  the  case  where  more  than  one 
element  is  required  to  fail  in  order  to  have  total  structural  failure.  This  method  is 
described  by  examining  a  four-element  structure  where  two  elements  must  fail  (see  Ap¬ 
pendix  C).  The  Pg.  values  are  determined  for  the  elements  by  one  of  the  methods  pre¬ 
viously  described.  The  resultant  Ps  for  the  four  elements  is  determined  from  applica¬ 

tion  of  the  series  approach.  It  should  be  noted  that  this  method  is  somewhat  less 
conservative  than  the  conventional  weakest  link  method.  It  is  possibly  more  realistic, 
especially  for  the  complex  states  that  exist  in  the  ATP  structure. 

NUMERICAL  RESULTS  AND  DISCUSSION 

Numerical  reliability  ( R )  results  are  listed  in  Table  2  for  the  fragmentation 
shell  subjected  to  an  internal  pressure  of  14  ksi  (proof-test  load)  as  a  function  of 

minimum  detectable  crack  size.  The  Warner  diagram  method,  as  outlined  in  the  text,  is 

represented  by  normal  stress-strength  PDF's  determined  from  FE  solution  and  laboratory 
strength  test  results.  The  FSS  results  are  obtained  by  applying  the  method  presented 
in  the  text  where  the  strength  PDF  is  the  result  of  using  the  simulation  scheme  (see 
Appendix  B).  The  Weibull  R  values  were  obtained  from  application  of  Equation  9.  It 
should  be  noted  that  both  Warner  and  Weibull  methods  do  not  show  variations  due  to 
changes  in  crack  sizes.  This  is  expected  since  they  do  not  consider  results  from  FSS 
procedures,  which  indicates  the  importance  of  finding  a  crack  of  at  least  0.025  inch  or 
less  to  establish  at  least  92%  probability  of  survival.  The  Warner  diagram  shows  an 
estimate  of  1  failure  in  1000,  and  the  Weibull  method  shows  an  estimate  of  222  failures 
in  1000.  A  coefficient  of  variation  of  10%  represents  the  variation  in  element  stress 
values  obtained  from  the  finite  element  solution. 

In  Table  3,  probability  of  survival  estimates  are  given  for  variation  in  crack 
sizes  of  0.020  inch.  The  effects  of  FE  solution  errors  (e.g.,  CV  values)  are  less  than 
for  larger  cracks  of  0.100  inch.  Mean  and  CV  values  for  Kjq  were  obtained  from  labora¬ 
tory  test  data.  The  material  yield  strength  of  140  ksi  provided  an  upper  bound  for 
calculated  ac  (allowable  stress)  obtained  from  the  FSS  approach. 

Table  4  provides  a  partial  listing  of  reliability  for  the  ATP  system  where  crack 
sizes  are  0.100  and  0.050  inch  with  variations  in  CV  of  5%  to  20%.  The  effects  of  FE 
errors  are  noted  as  in  Table  3:  small  variation  -  better  reliability;  large  variation  - 
poorer  reliability.  The  reliability  numbers  in  parentheses  are  the  results  from  appli¬ 
cation  of  the  series-parallel  method  described  in  the  text. 

The  series-parallel  method  which  requires  failure  of  all  adjacent  elements  to  have 
structural  failure  provides  a  much  less  conservative  estimate  of  reliability.  It  is 
possible  that  an  upper  and  lower  bound  on  reliability  of  this  structure  for  the  speci¬ 
fied  crack  sizes  could  be  a  series-parallel  system  and  the  weakest  link  approach, 
respectively.  With  tension  and  compressive  stresses  existing  in  this  structure,  it 
does  not  seem  advisable  to  consider  structural  failure  in  terms  of  only  one  given  ele¬ 
ment  failure.  It  is  also  unreasonable  to  assume  that  all  adjacent  elements  must  fail 
to  have  structural  failure,  particularly  if  a  bending  stress  exists  in  the  structure. 

Although  the  FSS  model  is  an  idealized  method  for  estimating  structural  reli¬ 
ability,  it  does  provide  a  desirable  alternative  to  the  present  fracture  mechanics 
approach  which  assumes  cracks  are  orientated  transverse  to  the  maximum  principal  stress, 
resulting  in  unnecessarily  high  rejection  rates.  The  ability  to  examine,  at  least 
qualitatively,  the  reliability  of  structures  as  related  to  the  ability  to  detect  flaws 
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or  cracks  of  various  magnitudes  can  provide  a  guide  for  future  NDT  development  pro¬ 
cedures.  If  more  information  was  known  regarding  structural  flaw  distributions,  the 
FSS  application  could  provide  an  excellent  reliability  tool,  certainly  one  superior  to 
the  present  laboratory  test  procedures  applied  to  brittle  materials.  It  should  be 
noted  that  in  laboratory  testing,  surface  flaws  are  often  removed  from  material  thereby 
preventing  an  accurate  representation  of  the  material's  strength  as  it  is  related  to 
the  structural  component. 


Table  2.  RELIABILITY  RESULTS  VERSUS  MINIMUM  DETECTABLE  CRACK  SIZE 


Minimum  Detectable 
Crack  Size  (in. ) 

Probability  of  Survival  - 

-  Shell 

(Proof-Test  Load) 

Warner  Diagram  Monte 

Carlo 

(F.M.) 

Weibull 

0.100 

0.999 

0.059 

0.778 

0.050 

0.999 

0.451 

0.778 

0.025 

0.999 

0.920 

0.778 

CV  (Element  Stress)  =  10* 
KjC  =  30  ksi/TnT,  CV  =  12% 
Yield  Strength  =  140  ksi 


Table  3.  STRUCTURAL 

RELIABILITY 

OF  SHELL 

(PROOF-TEST  CONDITION) 

Minimum  Detectable 
Crack  Size  (in.) 

CV  (Element  Stress) 

5% 

10% 

15% 

20% 

0.100 

0.210 

0.059 

0.027 

0.022 

0.075 

- 

- 

- 

- 

0.050 

0.557 

0.451 

0.378 

0.337 

0.020 

0.972 

0.966 

0.956 

0.935 

0.010 

- 

- 

- 

- 

Kjc  =  30  ksi /Tn. ,  CV  =  15% 


Yield  Strength  =  160  ksi 


Table  4.  STRUCTURAL  RELIABILITY  OF  ATP  SYSTEM 


CV  (Element  Stress) 
Crack  Size  (in. ) 

5% 

10% 

15% 

20% 

0.100 

0.353  (0.960) 

0.152  (0.841) 

0.079  (0.744) 

0.042  (0.569) 

0.075 

- 

- 

- 

- 

0.050 

0.606  (0.976) 

0.450  (0.970) 

0.378  (0.965) 

0.325  (0.954) 

0.025 

- 

- 

- 

- 

KIC  =  40  ksi/TnT,  CV  =  10% 

Yield  Strength  =  215  Tension,  Yield  Strength  =  240  Compression 


Numbers  in  parentheses  are  series-parallel  systems;  others,  weakest  link  concept. 
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In  application  of  the  FSS,  determination  of  the  proper  number  of  simulations  in 
the  Monte  Carlo  method  was  obtained  from  examining  the  convergence  rate  for  the  calcu¬ 
lated  reliability  numbers.  Instead  of  relying  on  some  elaborate  formulation  for  es¬ 
tablishing  the  proper  number  of  simulations,  a  chart  similar  to  the  one  shown  in 
Figure  8  was  used.  All  functional  parameters  were  increased  equally  in  number  to 
examine  all  effects  of  the  simulation.  To  examine  the  acceptability  of  this  method,  a 
comparison  was  made  for  R  using  the  Warner  diagram  approach  where  normal  PDF's  were 
calculated  from  a  prescribed  mean  and  standard  deviation  to  represent  the  stress-strength 
values.  Results  show  excellent  agreement  with  the  FSS  simulation  method,  using  the 
convergence  criteria  for  number  of  simulations  (see  Figure  8). 


Figure  8. 


CONCLUSIONS 

The  flaw  simulation  scheme  requires  that  each  structur 
fLaw  of  a  particular  geometry.  This  assumption  introduces 
limits  acceptability  of  the  quantitative  results  presented, 
exist,  it  has  been  demonstrated  that  application  of  this  sc 
analysis  in  relation  to  flaw  detection  capabilities  as  noted  in  Tables  2,  3,  and  4. 

The  Warner  diagram  and  Weibull  approach  provide  a  means  of  obtaining  structural 
reliability  without  application  of  the  conventional  fractural  mechanics  methodology. 
These  methods  can  provide  either  too  conservative  or  overly  optimistic  (e.g.,  Warner 
diagram)  reliability  results  as  demonstrated  in  Table  2. 

The  weakest  link  and  series-parallel  schemes  are  reasonable  computational  methods 
for  obtaining  reliability  error  bounds  as  related  to  structural  failure  mechanisms.  As 
noted  in  Table  4,  the  weakest  link  approach  represents  the  lower  bound  (conservative 
estimate)  while  series-parallel  determines  the  upper  bound.  The  latter  is  probably 
more  reaListic  for  the  ATP  system  with  the  other  more  applicable  for  shell  reliability 
computations . 

The  effects  of  increasing  uncertainties  in  application  of  the  FE  solution  are 
demonstrated  from  the  results  in  Tables  3  and  4.  The  results  in  both  tables  indicate  a 
reduction  in  reliability,  with  the  amount  of  decrease  depending  on  the  detectable  crack 
size . 

It  is  obvious  from  the  tabulated  results  that  in  order  to  have  acceptable  reli¬ 
ability  of  the  structures,  one  should  be  able  to  determine  relatively  small  size  cracks 
(e.g.,  0.025  inch  or  less). 


al  element  contain  an  isolated 
an  idealization  which  somewhat 
Although  these  limitations 
heme  can  provide  a  sensitivity 
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The  primary  contribution  of  the  flaw  simulation  scheme  is  the  ability  to  isolate 
the  effects  of  the  various  parameters  involved  in  determining  structural  reliability 
from  a  LEFM  technology.  The  ability  to  recognize  excessive  reduction  in  reliability 
due  to  relatively  small  errors  in  the  parameter  can  establish  need  for  improvement  in 
the  parametric  determination.  Accurate  quantitative  reliability  values  are  extremely 
difficult  to  predict;  therefore,  the  objective  should  be  to  reduce  the  error  in  the 
prediction  process  by  making  necessary  improvements  in  the  sensitive  parameter 
determinations . 

The  failure  rate  reported  by  the  fragmentation  shell  manufacturer  during  proof 
testing  was  from  2  to  5  percent  depending  on  the  production  lot.  Since  shells  con¬ 
taining  flaws  greater  than  0.020  in.  were  rejected  prior  to  proof  testing,  this  result 
correlates  exceptionally  well  with  the  results  shown  in  Table  3.  Although  the 
objective  of  the  paper  was  to  obtain  a  sensitivity  analysis  with  respect  to  the  critical 
stress  parameters,  the  results  have  indicated  possible  quantitative  acceptance  of  the 
FSS  reliability  computations. 
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APPENDIX  A.  PROBABILITY  OF  SURVIVAL  DETERMINATION 


s  =  Calculated  Stress 
S  *  Materials  Strength 

dPs  ■  (2(5!)  ds  /s  fy(S)  dS 

And  P$  ■  /dPs  -  J  tjlsif  J  fjISldslds 

Where  f(s!  and  f!S)  are  PDF  representation  for  stress 
and  strength  values  respectively. 


APPENDIX  B.  PROBABILITY  OF  ELEMENT  SURVIVAL 


M  M 

Ps«  *  1/M^  Z—i  a;  where  aj  *  I  l 

K  j,j  pj  1  '  (0  Otherwise | 

M  ■  Number  Simulations 

PsK  -  Probability  of  Survival  of  Element  K. 


APPENDIX  C.  SERIES-PARALLEL  CONCEPT  (TENSION-COMPRESSION) 

Assume  simply  four  element  structure 


H~?s,Tpsn^Tp^> 


Survival  Probability  for  Structure 


For  large  structures  all  elements  are  considered  with 
their  corresponding  adjacent  elements. 
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